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We report an inelastic-neutron-scattering and muon-spin-relaxation study of the effect of 2% spinless 
(Zn) impurities on the magnetic order and dynamics of YBa2Cu306,6j an underdoped high-temperature 
superconductor that exhibits a prominent spin-pseudogap in its normal state. Zn substitution induces 
static magnetic order at low temperatures and triggers a large-scale spectral-weight redistribution from 
the magnetic resonant mode at 38 meV into uniaxial, incommensurate spin excitations with energies well 
below the spin-pseudogap. These observations indicate a competition between incommensurate magnetic 
order and superconductivity close to a quantum critical point. Comparison to prior data on La2_jjSrjfCu04 
suggests that this behavior is universal for the layered copper oxides and analogous to impurity-induced 
magnetic order in one-dimensional quantum magnets. 



PACS numbers: 75.10.Kt, 74.72.Kf, 74.72.Gh, 74.25.Dw 

Spinless impurities have provided key insights into the 
spin correlations of low-dimensional quantum magnets. Re- 
search on quasi-one-dimensional (ID) quantum-disordered 
magnets with a spin gap has shown that every magnetic 
vacancy contributes to a Curie term in the uniform suscepti- 
bility and generates a cloud of antiferromagnetically corre- 
lated spins whose spatial extent is inversely proportional to 
the gap [ill . The striking discovery that a dilute concentra- 
tion of spinless impurities induces antiferromagnetic long- 
range order in a spin-Peierls system [2] can be attributed to 
constructive interference between spin-polarization clouds 
centered at different impurities. Related "order from dis- 
order" phenomena were later observed in a wide class of 
quasi-lD systems and shown to exhibit universal characteris- 
tics [3, 4]. In contrast to the comprehensive understanding 
that has been developed for ID magnets, theoretical work 
on spinless impurities in 2D quantum magnets has thus far 
been restricted to analytical calculations at quantum critical 
points [5-8] and numerical simulations for specific 2D spin 
Hamiltonians with quantum-disordered ground states [0 . 
Both approaches predict impurity-induced spin textures and 
order-from-disorder phenomena not unlike those observed 
in ID magnets. Experimental tests of these predictions have, 
however, proven difficult because of the paucity of physi- 
cal realizations of 2D spin liquids and the limited chemical 
solubility of nonmagnetic atoms in specific compounds. 

Layered copper oxides have been a fruitful testing ground 
for theories of 2D quantum magnetism and for the behavior 
of spinless impurities in two dimensions, despite additional 
complications introduced by the presence of itinerant charge 
carriers and d-wave superconductivity [fTl llOIL Spinless im- 
purities can be readily introduced into the CUO2 planes 
by replacing spin- 1/2 copper ions by nonmagnetic zinc or 
lithium. A large body of nuclear magnetic resonance (NMR) 
experiments on metallic and superconducting YBajCugOe+x 



(YBCOe+x) has demonstrated that such impurities generate 
a Curie term in the uniform susceptibility and a staggered 
spin polarization cloud on neighboring Cu sites, in close 
analogy to the ID systems discussed above [Jj . NMR ex- 
periments on other families of cuprates are much more 
limited, because intrinsic disorder due to dopant atoms 
broadens the NMR lines, so that the effect of Zn impurities 
is masked by native disorder llllll . Inelastic neutron scatter- 
ing (INS) experiments on both YBCOs.,.^ and La2_;cSr;^Cu04 
(LSCO) have revealed Zn-induced low-energy spin excita- 
tions, in qualitative agreement with the NMR results. In 
contrast to the ID systems [4], however, a universal picture 
of the impurity-induced magnetism has not yet emerged 
in the 2D cuprates, because the spatial character of the 
impurity-induced low-energy spin dynamics in both materi- 
als appeared to be rather different: commensurate antifer- 
romagnetic in YBCOg+j. I I121I13II and incommensurate (IC) 
in LSCO [14. 15]. Moreover, Zn-induced static incommen- 
surate magnetism, which was observed in LSCO II14H16II 
and interpreted as evidence of pinning of fluctuating stripes 
[]l7], has thus far not been detected in YBC05+^. For this 
reason, and because LSCO exhibits static IC magnetic order 
at some doping levels even in the absence of spinless impu- 
rities [18l |l9|] , the universality of this "order from disorder" 
phenomenon has been questionable. These considerations 
have motivated us to revisit the effect of spinless impurities 
on the magnetic correlations of YBCO5 6, a material with 
minimal native disorder that is well known for its large 
normal-state spin-pseudogap JT]] . 

The experiments were performed on an array of de- 
twinned, co-aligned YBa2(Cuo.98Zno.o2)306.6 (henceforth 
YBCOs etZn) single crystals of total mass 650 mg and su- 
perconducting Tj. = 30 K. The Zn content was determined 
by energy-dispersive x-ray analysis and inductively-coupled 
plasma spectroscopy. A twin domain population ratio of 



2 



4:1 enabled us to discriminate spin excitations propagat- 
ing along the a- and b-axes in the CUO2 planes, which is 
important in view of the recent observation of uniaxial, in- 
commensurate magnetism in YBCO6.45 [201121]. The INS 
measurements were carried out on the spectrometers INS 
(Institut Laue-Langevin, Grenoble), 4F2, IT , 2T (Orphee, 
Laboratoire Leon Brillouin, Saclay), and PUMA (FRM-II, 
Garching), using either a pyrolytic-graphite (PG) (002) or 
a Si (111) monochomator in combination with a PG (002) 
analyzer PG filters extinguished higher-order contamina- 
tions of the beam. The final wave vector was fixed to 4.1 
for excitation energies £ > 25 mey and to 2.662 or 
1.97 A~^ for £ < 25 meV Reference measurements were per- 
formed on a Zn-free YBCO5 5 array with = 61 K described 
previously [22] . Data on both arrays were compared by cal- 
ibration measurements on the optical phonon at 42.5 meV 
following previous work [23J, and by performing measure- 
ments under identical conditions. The wave vector transfers 
Q = [H,K,L) are given in units of the reciprocal-lattice vec- 
tors a*, b*, and c*. Q = (_H,K) stands for 2D projections of 
Q, and = (0.5, 0.5) is the wave vector characterizing 
antiferromagnetism in undoped YBCO5 q. 

Figure la,b shows INS scans taken along the a-axis of 
YBCOs.e^Zn at £ = 5 meV, well below the spin pseudogap of 
pristine YBCOe.e- A pair of well separated IC peaks symmet- 
rically displaced from is observed consistently in two 
Brillouin zones on top of an extraneous background arising 
from the sample mount. The background is independent 
of the sample orientation and can be readily subtracted 
(Fig. Ic). The signal is absent at L = and vanishingly 
small at large |Q|, and can thus be ascribed to acoustic spin 
fluctuations [23 J. Very similar peaks are also observed at 
£ = 7 and 10 meV (Fig. ld,e). As scans along the b-axis are 
peaked at Q^, the excitations are characteristic of uniaxial 
IC magnetism. Since low-energy spin excitations of com- 
parable intensity are not observed in control experiments 
on pristine YBCO5 6 (see below), they can be identified as a 
result of Zn substitution. 

The incommensurability of the impurity-induced spin ex- 
citations eluded detection in previous INS experiments on 
YBC05_,.^:Zn [12. .ISJ for three reasons: First, the signal 
from the previously investigated twinned crystals consists 
of equal contributions from two twin domains, so that the 
IC character is less apparent. Second, earlier studies were 
performed in the scattering plane spanned by the vectors 
Q = (1, 1,0) and (0,0, 1), which does not include the IC 
wave vector. Third, the effective momentum-space resolu- 
tion in the latter configuration is poorer than in the present 
experiment, because the long, vertical axis of the resolution 
ellipsoid is parallel to the CUO2 -planes. 

The INS profiles were fitted to Gaussians centered at 
Qj^p ± (5, 0) (lines in Figs. Ic-e). In the range £ = 5 — 10 
meV, both the incommensurability 5 = 0.10±0.05 and the 
intrinsic half-width-at-half-maximum AQ = 0.055 ± 0.025 
extracted from these fits were found to be energy indepen- 
dent within the error The incommensurability is consistent 
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FIG. 1: Constant-energy scans in YBCOg giZn at T = 10 K. (a,b) 
Uncorrected scans along (a) Q = (H, 0.5, 1.7) and (b) (H, 1.5, 1.7). 
The extraneous background was determined by rotating the sample 
by a) 30°, b) 10° and repeating the scans (open symbols), (c-f) 
Background-corrected scans along Q = (H, 1.5, 1.7) (full symbols) 
and (1.5, K, 1.7) (open symbols) at £ = 5, 7, 10, and 38 meV 



with the one that characterizes excitations above the spin- 
pseudogap of pristine YBCOg g ||22||. but significantly larger 
than the one observed in YBCO6.45 [120'] . This confirms that 
the Zn content and the doping level are independent control 
parameters [I j. Models according to which Zn impurities 
reduce the hole doping level [24>.25 J are inconsistent with 
our data. The intrinsic momentum width translates into 
a length scale of ^ = 2.9a, where a is the in-plane lattice 
spacing. This represents a lower bound on the intrinsic 
coherence length of the spin fluctuations, because inhomo- 
geneous broadening (arising from the doping dependence 
of 5 in combination with a possible slight nonuniformity of 
the doping level) may also contribute to AQ. In any case, 
the value of ^ determined in this way is in good agreement 
with NMR data [.261 and comparable to 5.8a, the mean 
distance between Zn impurities in the CUO2 planes. 

In view of the sharp excitation peaks, we have investi- 
gated the presence of magnetic order in YBCOg 5:Zn. Since a 
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FIG. 2: (color online) Zero-field muon-spin relaxation data for 
YBCOggiZn and reference samples YBa2Cu3055 and YBa2Cu306,45 
at r = 2 K. The lines are the results of fits described in the text. 



quasielastic neutron signal like the one previously observed 
in YBC06.35 (Ref. TJ) and YBCO6.45 (Ref. 'JT) could not 
be detected above background [28], we performed muon 
spin relaxation (/iSR) measurements, which are well suited 
to search for weak magnetic order. The experiments were 
performed with the GPS setup at the TimS beamline at 
the Swiss muon source S/iS at the Paul Scherrer Institute, 
Villigen, Switzerland. The time-resolved /xSR spectra of 
YBCOg giZn as well as reference spectra of Zn-free YBCO5 6 
and YBCO6.45 (Fig. 2) were fitted to the relaxation function 
P(t) = P(0)exp[(-At)''] X KT that consists of the product 
of a stretched exponential and a Kubo-Toyabe (JCT) func- 
tion, which account for the contribution of the electronic 
and nuclear magnetic moments, respectively. For YBCO5 6 
the relaxation is dominated by the nuclear moments, sug- 
gesting the absence of static or slowly fluctuating electronic 
moments. The YBC06.6:Zn spectra exhibit a more rapid de- 
polarization at temperatures below ~ 10 K (Fig. 2), which 
is a signature of dominant contributions from electronic 
moments. The absence of a precession frequency indicates 
that these electronic moments are either static and disor- 
dered or slowly fluctuating. The obtained values of ^ ~ 0.4 
for YBCO6.45 and § ~ 0.7 for YBCOg 6:Zn are characteristic 
of a spatially nonuniform distribution of the magnetic mo- 
ment amplitude or fluctuation rate II 1611 . These observations 
confirm prior /iSR results [JH |29l I30II according to which 
spinless impurities extend the range of magnetic order in the 
phase diagram of YBCOs+^f up to x ~ 0.7. The low-energy 
INS data (Fig. 1) and the close similarity between the /iSR 
spectra of YBC05 6:Zn and those of pure YBCO6.45 (Fig. 2) 
now indicate that this order is incommensurate, in contrast 
to the commensurate antiferromagnetism suspected on the 
basis of prior INS data on YBCOg+^iZn [12, 13]. 

We now describe the behavior of the spin system at higher 
temperatures and excitation energies. At T = 70 K (that is, 
above of both YBCOg.g and YBC06.6:Zn), low-energy IC 
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FIG. 3: (color online) Temperature dependence of the spin fluctua- 
tion intensity. (a,b) Scans of the dynamical magnetic susceptibility, 
j"(Q, E) (derived from the INS data by correcting for the detailed- 
balance factor), around Q = (0.5,0.5, 1.7) in Zn-substituted and 
Zn-free YBCO5 5 at energies and temperatures marked in the leg- 
end. The data were background-corrected as described in Fig. 1. 
The lines are the results of fits described in the text, (c) Ampli- 
tude of x" in YBCOg giZn at £ = 5 and 7 meV extracted from fits 
to constant-energy scans, (d) Amplitude of %" for YBCOg giZn 
(YBCOg.e) at £ = 38 (37) meV The lines are guides-to-the-eye. 



fluctuations persist in YBCO5 5:Zn without significant ther- 
mal broadening (Fig. 3a), and remain much more intense 
than those in YBCO5 6 (Fig. 3b). The onset temperature 
of the IC intensity in YBCOggtZn is ~ 150 K (Fig. 3c), 
comparable to the onset of the spin-pseudogap in YBCO5 5 
previously documented by NMR [ 1 ] . This shows that Zn 
substitution restores low-energy spin fluctuations not only 
in the superconducting state, but also in the pseudogap 
state. The characteristic temperature for the Zn-induced 
enhancement of the spin dynamics is also comparable to the 
onset temperature of IC spin fluctuations in YBCOg 45, which 
provides evidence for an electronic nematic state II20II . 

A converse spectral-weight shift is observed at higher en- 
ergies. Pure YBCOg 5 exhibits a sharp resonant mode with 
energy E^^^ = 38 meV in the superconducting state [22, 23], 
which is manifested by a sharp enhancement of the spin 
fluctuation intensity for E ~ E^^^ below (Fig. 3d) . While 
YBC05 6:Zn exhibits a similar narrowing of the magnetic 
response around and £ = 38 meV as pure YBCOe.e (Fig. 
10, followed by a dispersion away from at higher ener- 
gies (not shown), the superconductivity-induced anomaly is 
obliterated by Zn substitution (Fig. 3d) . This is consistent 
with previous INS studies close to optimal doping that had 
demonstrated progressive weakening of the resonant mode 
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FIG. 4: (color online) Energy dependence of the peak susceptibility 
x" SLtT = 2K extracted from constant-energy scans on YBCOg giZn 
and YBCOg 6- Normal-state data on pure YBCOg g are shown for 
comparison. The lines are guides-to-the-eye. 

with increasing Zn content without significant renormahza- 
tion of its characteristic energy [31 -33 J. 

Figure 4 provides a synopsis of the spin excitation spec- 
tra of YBCOg.e and YBC06 6:Zn. The YBC06 6:Zn spectrum 
at low temperatures is closely similar to the normal-state 
spectrum of YBCO5 5, except for the strongly enhanced in- 
tensity at low energies that goes along with the appearance 
of IC magnetic order. The low-temperature spectrum of 
YBCO5 6, on the other hand, is dominated by the sharp mag- 
netic resonant mode at 38 meV, which is a signature of the 
superconducting state [23j . Fingerprints of both forms of 
collective electronic order are thus clearly apparent in the 
spin excitation spectrum, and the spectral-weight redistribu- 
tion induced by Zn substitution demonstrates that spinless 
impurities nucleate static IC magnetism at the expense of 
superconductivity. This conclusion is consistent with the 
behavior expected in the vicinity of a quantum critical point 
separating both phases, and with prior /xSR studies that 
had indicated a local suppression of the superfluid density 
around Zn impurities [34, 35] . The two-phase competi- 
tion also explains the extreme sensitivity of the magnetic 
resonant mode to a minute concentration of Zn impurities 
II31H33II . which is difficult to explain in the framework of the 
random-phase approximation commonly used to describe 
the origin of this mode (13611 . 

In conclusion, our experiments on the prototypical spin- 
pseudogap material YBCO5 5 have shown that spinless im- 
purities induce incommensurate magnetic order and low- 
energy spin excitations with a correlation range comparable 
to the mean distance between the impurities. In conjunc- 
tion with prior data on LSCO 1 1141 115II , this indicates a 
universal order-from-disorder phenomenology for copper 
oxides with 2D electronic structure, closely analogous to the 
one previously established for quasi-lD quantum magnets 
|]2:-4j . In contrast to the latter systems, however, the spin 
excitations in the layered copper oxides are incommensu- 
rate. Our demonstration that the incommensurability in the 



YBC05+^ system is independent of impurity concentration 
and controlled by the density of itinerant charge carriers, 
as previously shown for LSCO [iHl |l5l |l8l [l9l , explains 
the difference to the Mott-insulating ID systems and fur- 
ther underscores the universality of impurity-induced spin 
correlations in low-dimensional quantum magnets. 
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